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Abstract 

A method is proposed for the experimental measurement of neutron separa- 
tion energies for nuclei far from stability. The procedure is based on deter- 
mining cross sections for the production of nuclei, by projectile fragmentation, 
for which only protons are removed but for which the number of neutrons is 
left unchanged. A simple Abrasion- Ablation analysis leads to a cross section 
prediction which is sensitive to the neutron separation energy after a single 
parameter is adjusted in comparison with data. Examples which illustrate 
the method are presented. 
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Few experimentally measured values of neutron separation energies have been determined 
for nuclei near the neutron drip-line [1, ?,?,?]. However, it is precisely in this region that such 
information is particularly necessary for investigating nucleosynthesis, and for testing models 
of nuclear structure at the greatest distance from the valley of stability [5,6,?]. As the values 
of TV and Z move away from the valley of stability on the neutron-rich side, the neutron 
separation-energy systematically is reduced [5,6,?] due primarily to the asymmetry in the 
proton-neutron composition. This asymmetry is reflected by the value of the parameter 5 
which is defined as (N—Z) / (N+Z). The separation energies for nuclei with large asymmetry 
are often extrapolated from information available close to the valley of stability [8]. This 
approach may be uncertain in describing how their values goes to zero with increasing S. 
However, this region tests the role of the symmetry energy most stringently in theoretical 
models, and hence provides information about the symmetry energy in the more general 
context of the nuclear equation-of-state [9] . Thus, the systematic behavior of the separation 
energy with increasing values of 5 offers direct evidence for such effects. 

One of the most successful methods for the production of neutron-rich rare nuclides has 
been the fast projectile fragmentation process [10,11]. This method has been used to extend 
the list of particle stable nuclei to the extremes [12,?]. In most of the cases where the values 
of the neutron separation-energy have been measured, the values are in excess of 5 MeV. 
In the case of nuclides very close to the drip-line the pairing energy may cause 2-neutron 
separation energies to be less than the 1-neutron separation energies. For those cases the 
lowest separation energy, referred to as S in this article, be it for one or two neutrons, is 
the one of interest. In this paper we suggest that, under certain circumstances, the same 
measurement which provides the verification for the existence of a rare nucleus may also be 
used to estimate S. 

One avenue to this information lies in the recent suggestion that "cold" fast fragmentation 
[10,11] seems an efficient method of producing these extremely rare nuclei. The simple sce- 
nario for this process follows the concepts of the Abrasion-Ablation (A-A) models [14,10,11]. 
From that point of view, a direct reaction (abrasion) removes a number of nucleons, leaving 
the residue excited, and free to lose more nucleons by evaporation (ablation). 

To minimize the uncertainties associated with the abrasion and ablation processes, we 
focus on the production of fragments where only protons and no neutrons are removed from 
the projectile, i.e., the abrasion process removes the protons, and leaves the residue with 
too little excitation energy to permit further loss of particles (neutrons). Such a production 
mechanism is referred as a "p-removal chain" in this article. If the residue is neutron-rich, 
it decays by neutron emission. Thus the upper limit on the excitation energy is the neutron 
separation-energy, S. It is this feature that provides the production cross section with the 
sensitivity to the separation energy. By limiting our attention to the nuclei produced in 
p-removal chains, we avoid the ambiguities related to specific evaporation models. 

Using the frame work of the A-A model [11], the cross section to produce a nucleus with 
(Z — x) protons and N neutrons from the fragmentation of a projectile with Z proton and N 
neutrons can be written as a x = Abr x ■ Abl x . Here, the factor Abr x gives the cross section for 
removing x protons (and no neutrons) by abrasion, and Abl x is the dimensionless probability 
that the residue will not further decay following the removal of those x protons. 

The factor Abr x can be estimated by the geometric overlap [15] of projectile with the 
target. This can provide the cross section for the removal of x particles [14,10]. This 
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cross section must then be multiplied by the probability that all the abraded particles are 
protons. By assuming that the positions of neutrons and protons in the projectile are 
uncorrelated, the probability that all x of the abraded particles are protons can be written 
as (Z\/(Z — x)\)/((N + Z)\/(N + Z — x)\). Both of these assumptions (geometric overlap and 
uncorrelated positions) are simplistic but their validity can be calibrated by comparison with 
a measured set of cross sections for projectile fragmentation where the separation energies 
of the resulting nuclei are known. 

The crucial factor is Abl x which depends on both the distribution of excitation energy 
following the removal of x protons, F X (E*), and the separation energy, S x , for the nucleus 
produced by this removal. Specifically, Abl x is the integral of the excitation function from 
zero to the separation-energy. 

Clearly the form of the excitation function is a critical component of this procedure. 
The literature of A-A models suggests that this distribution function may be quite uncer- 
tain [14,10]. Recent approaches [10,?] suggest that the distribution function, F X (E*), is a 
convolution of x distribution functions, /i(e*), where /i(e*) is the function for the removal 
of a single nucleon: 

WE*) = I Il(de* /i(e*)) 5(±e* - E*) (1) 

J i=l i=l 

The functional form of /i(e*) is, however, not well determined, and there is great uncertainty 
as to the mean value of the excitation energy, < e* >, it provides. Some of this uncertainty 
can be removed by fitting calculated cross sections to sets of measured cross sections. The 
fitting would be accomplished by the adjustment of < e* >. 

One form of the suggested single particle excitation distribution, /i(e*), widely used [14] 
is the "triangle" distribution. This has the form /i(e*) = 2/E m (l — e*/E m ) for e* < E m 
and the average excitation energy, < e* > is E m /3. A wide range values of E m have been 
suggested in different models [14,?]. Convolution of this (triangle) single-particle distribution 
leads to a value for Abl x which is approximately (2S x /(3 < e* >)) x /x\, for S x « 3 < e* >. 
Exact values for Abl x can be obtained with 

Abl x = C tri (x) ■ (2S x /(3 < e* >)) x /x\, (2) 

where 

Ctn(x) = E(-S*/(3 < e* >)) s ■ (x\ 2 /(s\(x + s)\(x - s)\)) (3) 

The value of C tri (x) goes to 1.0 for small values of (S x /(3 < e* >)), and Abl x is seen to be 
a function of the parameter (2S x /(3 < e* >)). 

We have also considered a different form for the single-particle excitation function, i.e., 
the exponential function /i(e*) = 1/ < e* > exp(— e*/ < e* >) with an average excitation 
energy < e* >. A convolution of this function provides an x-particle distribution function 
of the form 

F X (E*) = (E*/ < e* >) x ~ l /(x - l)!exp(-E7 < e* >) (4) 
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When this function is integrated from zero to the separation energy S x , one obtains a value 
for Abl x which is approximately (S x / < e* >) x /x\, and an exact expression can be calculated 
as a function of (S x / < e* >): 

Abl x = C exp (x) ■ (S x / < e* >) x /x\, (5) 

with 

C exp (x) = £(-^/ < e * >) s (^/((^ + s ) s 0) (6) 

The value of C exp (x) goes to 1.0 for small values of (S x / < e* >). 

For small values of S x , the functional form of the Abl x (as a function of different pa- 
rameters) is the same for the two single-particle distribution functions. For the "triangle" 
distribution the parameter is (2S x /(3 < e* >)), while for the exponential distribution the 
parameter is (S x / < e* >). Hence one might expect that a choice of < e* > in the "triangle" 
distribution which is 2/3 the choice of< e* > in the exponential distribution would predict 
similar cross sections. When the full x-particle excitation distributions are used, however, 
this scaling is only approximate. 

For the nuclei in a given p-removal chain, estimates of unknown values of S x can be 
extracted from measured values of cross sections, a x , after a value of < e* > has been 
adjusted to fit the data for nuclei with known separation energies. 

In illustrative examples below we have found that we can, generally, well represent the 
data by first calculating Abr x using simple assumptions, and then adjusting the parameter 
< e* > to give Abl x . The assumptions for Abr x include both the estimation of the removal 
cross section from the geometric-overlap of target and projectile, and also the calculation 
of probability for obtaining pure proton removal by assuming uncorrelated positions for the 
nucleons. Any systematic correction required in Abr x may possibly be accommodated by 
adjusting < e* > in Abl x . The ambiguity in the choice of the single particle distribution 
function, f\, prevents a unique determination of the mean excitation energy through the 
fitting. However, we do find interesting systematic changes in the required values of < e* > 
which appear to depend on the mass and isospin values of the fragmenting projectile. 

We have considered, for illustration, the p-removal chains given in reference [10] for 208 Pb. 
197 Au and 136 Xe. We have also examined data in the literature for m Kr [16], and 48 Ca [12], 
and, in addition, preliminary results for 58 Ni which is under current investigation [17]. We 
first tested the form for the cross section suggested in the expressions of Eqs. 2 and 5, by 
fitting to data for the fragmentation production of nuclei where the separation energies are 
known. Both the "triangle", and the exponential forms for the excitation functions were 
used, and fits to the data were achieved by adjustment of the respective values of < e* >. 

In Fig. 1, we plot the fitted cross sections for the fragmentation of 86 Kr. The separation 
energies [4] are known for all first 5 members the p-removal chain (the data only covers 2-5). 
There are 3 degrees of freedom for this fit providing a \ 2 P er degree of freedom of .99 for 
the exponential distribution. The respective fitting values of < e* > are 11.7 MeV for the 
"triangle" distribution, and 16.6 MeV for the exponential distribution. 

In Table I we list the values of the x 2 / degree of freedom for fits to other data sets including 
only nuclei where the separation energies are known [4]. Excellent fits can be achieved. 
We have also listed values of the fitting parameters < e* > with estimates of deviations 
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(both plus and minus) for 70% confidence. There are variations in the values of these 
fitting parameters from one reaction to the another. In each case an approximate ratio of 
approximately 2:3 is found for the values related to the "triangle" and the exponential /i(e*) 
distributions with practically no difference in the resulting fits from the two distributions. 
Except for the Xe fragmentation, the data and calculated values, based on the excitation 
energies which provide the best fit [18], are plotted in Figures 1 to 4. Preliminary results 
show similar behavior for the p-removal chain of 58 Ni for which the chain has been measured 
through 8 p-removals [17]. 

While the current data is quite limited we also examined the predictive power of the 
method for nuclei where the separation energy is unknown, even lacking estimation by 
extrapolation. For this we looked at 204 Pt which was measured in the chain from 208 Pb 
projectile [19]. In Fig.3 we show the best fit to the data for the first three nuclei in the 
chain, where separation energies are known or estimated [4]. In the inserted graph we show 
the sensitivity to the assumed binding energies of the 204 Pt using the values of < e* > 
which best fit the first three members of the chain given in Table I. The apparent estimate 
for the separation energy is about 5 MeV with large uncertainty due to the experimental 
uncertainty in the measured cross section for 204 Pt. This value is well in line with systematic 
decrease of the separation energy with increased (N — Z). 

We have also attempted to estimate the separation energy of 41 Al which has recently 
been observed [12] for the first time. Even though the fragmentation of neutron-rich 48 Ca 
[12,13,20] has been studied intensely in the past few years, there is no systematic measure- 
ment of cross sections for the p-removal chain. The separation energies are not known for the 
nuclei with more than four protons removed. Some of the extrapolated values have uncer- 
tainties of much more than 1 MeV [4,?]. Thus in principle the 48 Ca would be a good place to 
fully examine our method. We have examined the existing data to make a rough estimate. 
We found the fragmentation cross sections have been measured from previous studies [20] for 
two nuclei, 45 CZ and 44 Si, which have 3 and 4 protons removed from A8 Ca. This experiment 
used a target of 9 Be. Using the separation energies 6.241 MeV and 5.21 MeV [4,?] respec- 
tively for 45 Cl and 44 Si the fit parameters listed in Table I are obtained. We next examined 
the results from the experiment [12] which first observed the 41 Al nucleus corresponding to 
the removal of 7 protons from the projectile. Unfortunately, these data were obtained with 
a 181 Ta target. To connect this point with the other two points (obtained with a Be target), 
we used abrasion calculations which suggest that the difference in targets provides a cross 
section from Be which is 0.545 times the value obtained with a Ta target. The reduction 
arises primarily from the difference in the respective size of the impact parameters for the 
two targets. We plot in Fig. 4 a point for x = 7 ( 41 Al) at a cross section of 4.4 pb, which 
is the value of 8 pb [12,?], reported for the Ta target, scaled down by the estimated ratio 
of cross sections. (We have also scaled down the error bar.) The insert in Fig. 4 shows the 
dependence of the calculated cross section for x = 7 as a function of the separation energy. 
The experimental uncertainty is high since only three events were observed [12,?]. Even 
so, extraction of a value of 3.5 ± 0.5 MeV would be consistent with the information in the 
inserted graph. We can not claim that this value is, indeed, the separation energy of 41 Al 
due to the fact that two different targets were used, and there is a scarcity of information 
in the p-removal chain (For example, there are no measured cross-sections and no accurate 
separation energies for for x=5 and 6 ( 43 P and 42 Si) isotopes). However, this exercise shows 
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the potential for extracting the separation energies for AA S, A3 P, 42 Si, as well as 41 Al, the 
nuclei with 4, 5, 6, and 7 protons removed from 48 Ca. This might be accomplished with 
careful measurements of the complete p- removal chain from x=l to 7 with one target and 
one beam energy. Specifically, additional data in the x=l-3 region where the S x values are 
known by observation will provide greater constraints on the values of < e* >. 

Finally, we have examined the situation for the fragmentation of 197 Au where a chain of 
5 proton removal is reported in ref. [10] for a target of 9 Be (solid points in Fig. 5. For the 
fragmentation of the 197 Au using an 27 Al target nuclei, there are p-removal cross sections 
up to x=3 [10,22]. For the case of 197 Au projectiles with 27 Al and 9 Be targets, abrasion 
estimates suggest a 5% reduction in going from the larger to the smaller target [23]. The 
three open points in Fig. 5 are the 27 Al data scaled down by 5%. They are consistently 
higher than the corresponding 9 Be data (solid points). If we apply the fitting procedures to 
this set of data, using the values of < e* > listed in Table I, we obtain a separation energy 
greater than 7 MeV for both 193 Re, and 192 W nuclei. These values are clearly inconsistent 
with systematic trends and expectations, both of which would have led to values below 7.0 
Mev. For comparisions, the solid and dashed lines are calculations using the best fit < e* > 
listed in Table I for 9 Be (36.3 MeV) and for 27 Al (32.2 MeV) targets, with the assumption of 
an exponential energy distribution. The upper and lower curve in each pair of lines use the 
separation energy of 7.0 and 6.5 MeV respectively as the separation energy for both 193 Re 
and 192 W nuclei. The calculated cross-sections are lower than the experimental values. In 
brief, the reported cross sections for the 197 Au + 9 Be reaction do not lead to reasonable 
separation energies for the last two members of the proton chain. The reasons for this 
failure are not clear at this time. 

In summary, the illustrated calculations show that excellent agreement with fragment 
cross-sections, where the fragment neutron separation energies are known, can be obtained 
when the simple estimates are used with the Abrasion-Ablation model, and a single param- 
eter, < e* >, is adjusted. The quality of agreement is equally good for both the "triangle" 
and the exponential single particle excitation distribution functions. The two distributions 
require values of the mean energy, which are approximately in the ratio of 2:3. The smaller 
the parameter the slower the fall of cross section with the number of protons removed. The 
quality of the fit inspires confidence in the use of the A-A model for calculating the p- 
removal chains. Once the parameter is determined for each chain the only remaining input 
is the set of separation energies. From some of the data in the literature we were able to 
suggest the power of the p-removal method for observing unknown separation energies in 
204 Pt and 41 Al. A puzzling disagreement was found for the unknown separation energies of 
193 Re and 192 W in the chain reported for the fragmentation of 197 Au when the target was 
9 Be. Clearly, more data and more understanding of the uncertainties in the cross-section 
measurements are needed to confirm the utility of the method. Since the procedure may 
be generally applied to all p-removal chains, it opens an avenue for measuring separation 
energies for neutron rich nuclei near the drip-line as illustrated by the fragmentation of 
extremely neutron-rich projectiles such as A8 Ca. While the method can not compete with 
dedicated mass measurements where masses can be measured to uncertainties better than 
10~ 7 [5-7], the simplicity of cross section measurements with fragment separators may allow 
the wide use of this method to measure the separation energies for extremely neutron rich 
nuclei to a couple hundred keV as this energy decreases toward zero at the drip-line. 
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TABLES 



TABLE I. Values of < e* > are obtained by fitting the measured cross sections of nuclides with 
know separation energies. The columns labelled + and - indicated the deviations of the the best 
fit values with 70% confidence. The 2 nd and Q th column reflect the goodness of the the fit for the 
"triangle" (tri) and exponential (exp) distributions. 



Reaction 


X 2 /dof 
tri. 


< e* > 
tri. 




+ 


X 2 /dof 
exp. 


< e* > 
exp. 




+ 


208 Pb + Cu [19] 


0.38 


18.4 


1.1 


1.5 


0.42 


26.6 


1.8 


2.1 


197 Au + 27 Al [22] 


0.87 


22.4 


1.6 


3.6 


0.88 


32.2 


3.8 


5.2 


197 Au + 9 Be [10] 


1.87 


25.0 


1.4 


1.8 


1.58 


36.3 


2.2 


2.6 


136 Ie + 9 Be [22] 


0.36 


23.8 


2.8 


2.6 


0.36 


34.2 


3.8 


5.6 


86 Kr + 9 Be [16] 


1.45 


11.7 


0.3 


0.25 


0.99 


16.6 


0.4 


0.45 


48 Ca + 9 Be [20] 


1.24 


7.70 


0.35 


0.4 


1.81 


10.80 


0.45 


0.60 
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FIGURE CAPTIONS: 

Figure 1. Measured cross sections (solid circles) for the production of p-removal nuclides 
with N = 50 , m ~ x Z from the fragmentation of 86 Kr with a target of 9 Be [16]. Lines are 
predictions described in text using two different excitation distributions with the adjustment 
of a single parameter, given in Table I. 

Figure 2. Measured cross sections for the production of p-removal nuclides with N = 
126, 208-1 Z from fragmentation of of 208 Pb with a target of m Cu [19]. Lines are predictions 
described in text using two different excitation distributions with the adjustment of a single 
parameter, given in Table I. The insert shows the predicted cross-sections as a function of 
the separation energy for 204 Pt nuclei. The horizontal solid and dashed lines are measured 
cross-sections. 

Figure 3. Measured Cross sections for the production of p-removal nuclides with 
N = 28, A8 ~ X Z from fragmenation of A8 Ca with 9 Be target [20] (solid points). Lines are 
predictions for x=l,2,3,4 described in text using two different excitation distributions with 
the adjustment of a single parameter, given in Table I. The open point for x=7 ( 41 Al) is 
obtained from separate experiment with 181 Ta target [12] and adjusted as describe in the 
text. The insert shows the predicted cross-sections as a function of the separation energy 
for 41 Al nuclei. The horizontal solid and dashed lines are measured cross-sections. 

Figure 4. Measured cross sections for the production of p-removal nuclides with N = 
118, 197 ~ X Z from a projectile of 197 Au with 9 Be target [10] (solid points). The open circle 
are data [22] for 197 Au + 27 Al scaled down by 0.95. Lines are predictions described in the 
text. 
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Cross Section (b) 




